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1. INTRODUCTION 


The model Hamiltonian most frequently used in literature 
for description of the chemisorption phenomenon is the 
Hamiltonian introduced by Anderson [1] and used by Newns [2) 
and Edwards. and Newns {3) - so called Newns-Anderson_ (NA) 
Hamiltonian. The Anderson model has first been introduced to 
describe a magnetic impurity in a metal but also is very useful 
to formulate fundamental microscopic foundations of tne 
chemisorption theory. However, the NA model does not take into 
account many 'important effects which may lead to considerable 
qualitative differences in some chemisorption characteristics 
(41. The model considered by us in this paper ıs a direct 
generalization of the Newns-Anderson mode: to the case of final 


318 E. Taranko, R. Taranko... 


coverages of randomly arranged adatoms [5]. The problem of 
constructing more realistic model of chemisorption has 
thorouglily been studied by us for the case of one admixture in 
Ref. [6). At present, the microscopic theories of the 
chemisorption systems can be roughly divided into three classes 
depending on methods used for the calculations of various 
chemisorption characteristics.In the first class we can place 
the methods using the local density functional method (LDF). 
Apparently, the description of the coverage dependence of the 
properties of the chemisorbed overlayers was for the first time 
realized within LDF {7}. In the framework of this scheme the 
adsorbate ion cores were modelled by a thin jellium slab 
adsorbed on a jellium surface .iIn the next step there 
appeared extensions of this approach, see, for example Refs. 
[8. 9], where an approximate description of the metallic 
substuste DY Set Der ASEM nE ne ee T SYS CEN or MY 
including the substrate structure within f2MSULerde:; 
perturbation theory for jellium plus adatom system was used. In 
Refs el LOT) Ehe alkali-metal adsorption on metal surfaces 
(the alkali-metal-overlayer-jellium system ) was investigated 
and a fùúliy self consistent electronic structure was 
obtained. In Ref. [14] a further improvement was achieved, i.e. 
à coverage dependence of the electronic structure of 
alkali-metal adatoms on a metal surface was studied treating 
both the adatom and substrate as discrete atoms. Although the 
"ab initio" method give very useful information, it involves a 
large amount of computer calculations and all the calculations 
must be repeated for new sets of an input data (for example, 
for a value of the coverage rate). In addition, these 
calculations rather reproduce the experimentally observed 
features and not explain it. The second class of methods 
consists essentially of .the quantum chemistry approaches (see, 
e.g. cluster methods {15]) and the third class comprises the 
model Hamiltonian methods (see e.g. 1[4]). Especially, the 
method of model Hamiltonians can be well suited for better 
understanding and not only for reproducing of the experimental 
results.It is to be noted that model approach is wide open to 
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criticism. However, in comparison with the "ab initio” 
calculations it is an indispensable tool for explanation of 
basic mechanism of physical processes, which essentially 
Promotes their understanding. 

In the present paper we consider the problem of many 
adatoms simultaneously existing on the metal surface within the 
model Hamiltonian method. The model Hamiltonian used here was 
recently obtained by us generalizes the Newns-Anderson 
Hamiltonian. The electron correlation effects are included 
within the second-order self-energy matrix formalism of Brenig 
and Schónhammer [16]. We also discuss here general solutions 
obtained within the equation of motion method for a case of 
submonolayer coverages. 

The paper is organized as follows. In the section 2 we 
give the approximate model Hamiltonian intended to describe the 
chemisorbed hydrogen overlayer on a model transition metal 
surface and compare this model with others known from the 
literature. Im the section 3 we present shortly some methods 
which incorporate the electron-electron interaction (on the 


adsorbate atom) and generalise the Hartree-Fock solutions. 
2. THE BASIC MODEL HAMILTONIAN.ONE ADATOM CASE 


Recently, we have obtained from the microscopic 
considerations the Hamiltonian intended to describe the system 
of adatoms located on the metal surface [5]. The system adatom 
plus metal surface is described by the Hamiltonian (we 


consider the hydrogen adatom in the 1 s ground state): 


hal L2 N+ Nel NH 
s3 am tE 2 VE) - > V(R, AIR PRU. (n 


Here the potential V(X-Y) describes the Coulomb interaction 
between electrons placed at the-positions X and Y, P(X) is the 
electrostatic interaction of electrons with t' electric fields 
of the metallic substrate. The constant L3 repre sents the 


energy of electrostatic interaction between ions of the metal 
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substrate and between these ions and the proton of the adatom. 
After performing the second quantization procedure (for details 
see Ref. (5]) we obtain: 


= €. > —--— wt — 
H 2. Eanna Un, n, st oe inget k kite MT. ka Tag Oeo Oe 
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elements calculated for the Coulomb interaction V(k-¥) between 
states described by the substrate electron wave functions a é 
pp. p? E $? and adatom electron function Pa and $? ; oy D p? ; 
adat TUR Other abbreviations have their OAM Ko 
(see, also [5]). Very similar model Hamiltonian, although 
without some terms present in (2) was obtained by Doyen [17] 
and Doyen and Ert! [18, 19]. Starting with the  Newns-Anderson 
Hamiltonian and adding additional self-consistent fields 
arising from electron-electron interactions they obtained 
parametrized model Hamiltonian which treats the interaction 
between adsorbate and metal electrons in some details. In our 
case we have started from very general physical description - 


Eq.(1), and as a result we have obtained the Newns-Anderson 
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Hamiltonian plus many other terms describing in detail the 
chemisorption process. To some extent, we can give a very rough 
physical interpretation for terms in the Hamiltonian - Eq. (2) 
(for comparison, see. Ref. [17]). For example, the fourth term 
describes repulsion between electrons localized on adatom and 
substrate electrons, the next term corresponds to the 
attraction of metal electrons by the positive ion core of the 
adatom and so on. Very important for our further considerations 
is the seventh term which can be interpreted as an interaction 
of overlap charge; with electrons on adatom.In other words, this 
term describes the influence of the fractional occupation of 
the surface impurity resonance on the electron transfer between 
the adatom and substrate metal. It is also important to note 
that for a consistent treatment of the polarization properties 
of the metal substrate it is insufficient to characterize its 
electronic properties by the parameter e (the energy spectrum 
of the substrate electron band), but it is necessary to 
introduce a long range part of the Coulomb interaction between 


substrate electrons v Be mno {20], or introduce plasmons [21]. 
Corresponding terrm is escent in our Hamiltonian. 

In the following we confine ourselves to the model in 
which together with the standard interactions present in the NA 
model we include only effects connected with the influence of 
the adatom orbital occupation on the charge transfer between an 
adatom and a metal substrate. In this case the Hamiltonian 


reads as [5, 6]: 


Ak 


t6 T i 
+ 
7 Vana Ta- ) Vag Gee * h.c.] ; 


The last term can be simplified as for 18,1 » 0.5 À we obtain 
MM. ^ aV . where the parameter a depends on the distance 
between the adatom and the metal surface (0S aS 1). We have 
checked in Ref.[6] that this approximation works relatively 


good. In Ref. [5] we have performed the self-consistent 
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calculations of the adatom electron charge and chemisorption 
energy in the framework of the Hartree-Fock approximation, i.e. 


for Hamiltonian 


inte 2260" te * 27 54 FLW aan -« SV RECI EIPS 
EA NET -«€n,. Pye Rp Qge* h.c] - (4) 


-U£n, >< nar ee (Vig C4 o Cg ge? * h.c.). 


The calculations were performed for parameters modelling the 
hydrogen chemisorption on transition metal surfaces. The 
results indicate that the charge transfer to the adatom is 
considerably damped with increasing value of the parameter a. 


This means that the additional term -V aa? n, ,a,jàgy, in NA 
Hamiltonian, which 1n a Hartree-Fock approximation leads to 
spin dependent effective hopping matrix element and 
renormalizes in special way the energy of an electron moving in 
the field of the adparticle ion core, plays an essential role 
in description of the adsorption system. Note, that the 
semiempirical one-particle Hamiltonian obtained for description 
of the chemisorption process in Ref.(19] leads to some 
difficulties . In this paper Doyen and Ertl checked that 
reasonable agreement with experiment could be only achieved if 
the hopping matrix element for different spins were allowed to 
differ. In this case it becomes obvious that the corresponding 
Hamiltonian is not the Hartree-Fock operator of a proper 
many-particle Hamiltonian. In order to.improve this situation 
one has to add to the Hamiltonian the term like ne Sar ca nre 
+ h.c., which was absent in model explored in Ref.[19]. In 
our formulation this term is present in the Hartree-Fock 
Hamiltonian from very beginning. It is contained in the second 
term in Eq.(4) leading to additional spin-dependence of the 
effective adatom energy level position (through the correlation 


function <a* a» >). 
A? ko 
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In the analysis we have performed above, the 
electron-electron interaction was treated in the framework of 
the Hartree-Fock approximation. It is difficult to take into 
account the correlation effects because the small parameter 
does not exists in the theory. This means that one has to sum 
an infinite series of the most divergent diagrams or construct 
the Green Function with the approximate self-energy operator. 
This solutions should be. of course, identical with the known 
ones for some special case. A great number of papers have been 
devoted to the study of the correlation effects in the 
chemisorption theory beyond the Hartree-Fock approximation. 
Here we mention merely few of them. First of all we have 
works by Brenig and Schónhammer [16], Yosida and Yamada [22], 
Anda, Majlis and Grempel [23], Sebastian and Rangarajan [24], 
Bell and Madhukar [25], Lacroix [26] and by Munoz et al. (27). 
Especially two of them are interesting. First, it is work by 
Brenig and Schónhammer [16] in which the correlation effects on 
adatom are treated within so called second order (in V - the 
single particle hopping strength between the adatom and the 
substrate atom) self-energy matrix formalism. This method was 
also transferred with success by Schónhammer (28] to 
considerations of the electron correlation effects in the 
standard many-body Hubbard Hamiltonian. Second work is the one 
by Munoz et al. [27]. In this case the self-energy expression 
interpolates between two opposite cases, i.e. between the 
second-order perturbation formula and formula for strong 
correlation limit. Here we use the second-order self-energy 
matrix formalism described in [16] to consider the electron 
correlation effects within our generalized NA Hamiltonian, i.e. 
with the term 


25 Nu Nyc lag Aget h.c, (5) 
kG 


included. Using the double-time retarded Green functions [29] 


builded up from the operators n, ,a,, anda ,(1- n, ,) 


we can try to solve the Dyson equation for this Green function 
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up to second order in Powers of Vp (here we used VA APA aV È 
[6]).As a result for the self-energy expression we Snai [30] 


Na-g? Lnag? 0 zi 
ver ena jun en ea) 


4- T,-e? O 4-0, e) (6) 
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and f(E) is a Fermi function. 

In Ref. [30] we have performed the self-consistent 
numerical calculations of the adatom charge for model hydrogen 
on titanium surface: chemisorption system for different values 
of the parameter a. The parameters for calculations were taken 
from Ref. [2]. The effect obtained after introduction of the 
self-energy corrections can be described as follows. For small 
values of a we have much smaller values of the adatom 
electronic charge comparing with the results obtained for 
standard NA model. For moderate values of a, i.e. fora ^ 0.5 
(we remember that a and distance between the adatom and 
substrate's surface depend on each other) the results obtained 
in Hartree-Fock approximation, as well as beyond Tipula 
approximation are comparable. For greater values of a the 
charge transfer "stabilizes" around the neutral adatom case 
more and more positive charged. 

On the other hand, we expect on physical grounds, that for 
a large separation distance (for large a) the hydrogen adatom 
shouid be in a nearly neutral state. Thus we can say, that the 
generalized NA Hamiltonian (i.e.,NA Hamiltonian plus correction 
expressed by Eq. (5)) better describes this chemisorption 
system. This description is sufficiently good for moderable 
values of a even in the Hartree-Fock approximation. From the 
comparison with results obtained for whole range of the 
parameter a it appears that generalized NA model combined with 
the self-energy corrections obtained within the method of 
Brenig and Schonhammer can be helpful in understanding of the 
chemisorption process. 


3. CHEMISORPTION AT FINITE SUBMONOLAYER COVERAGES 
The extension of the one-adatom chemisorption case to the 


submonolayer coverage can be described in analogous way as 
follows [31]: 
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The Hamiltonian in the form presented here contains the 
contributions which play the most important role in describing 
the electronic properties of adlayers. The other terms not 
written explicitly here are contained in R and involve 
processes which are not usually considered in simple 
chemisorption models. These terms describe the interaction of 
the valence adatom electrons with the redistributed charge 
cloud in the substrate metal induced by the adatom electrons 
themselves, effects describing the correlations between adatom 
charge fluctuations and the redistributed substrate charge 
density and others effects of. much higher order (for details see 
Refs {5, 311). Having in mind such complicated system as a 
chemisorbed laver on the metal surface we restrict ourselves in 
the first step to the simplest case,i.e. we consider the model 
described by the Hamiltonian: 


H= a E; ET € Naast S » Na Tag Ne -¢ * 
s (8) 


+> Veg NO se Qge* h.c. 
ka 


Model (8) is a generalization of the well-known NA Hamiltonian 
to the case of a stochastic arrangement of adatoms on the metal 
surface. The sum over o is carried out over all the adsorption 
centres. The operator Ny = ce has eigenvalues equal to 0 -or 
1, and the operators Ge and Ca are the Fermi amplitudes of the 
creation and annihilation of an adsorbed ion (adion) at the 
adsorption center a. We consider only the case when there is no 


more than one adatom ion in each adsorption center. At complete 


Electron Correlation Effects. . . 


monolayer coverage the adatoms form a two-dimensional lattice 
with the periodicity commensurate with that of the substrate 
surface. For the coverage less than unity the two-dimensional 
chemisorbed layer is considered as a substitutional alloy 
composed of two types of species, i.e. adatoms and vacancies. 
In the following we are going to consider the 
electron-electron correlation effects within the method of 
Brenig and Schonhammer (16] and the adatoms properties related 
with its position on the substrate surface are involved in 
correlators of the Ising type «N.N n, ,n, .>. At this stage of 
our calculations, in order to introduce | simplifications in 
handling higher-order Green functions composed of ion and 
electron operators we can use decoupling approximations. One of 


them concerns the Ising type correlators and here we proceed 
Las follows 


{Na Np? = bap Na? * (17855) ON € No) : (9) 


mo 


“1 1 4 4 
ENERGY ° ENERGY 

Fig.1. The adatom local density of states for different 
coverages (first curve in each panel corresponds 
to @=0, second for 60:5 and last for 6=1.0 ) in 
the symmetric case for U/W=0.15 (left panel) and 
U/W=0.3 (right panel ) . Energy in units of the 
substrate half bandwidth and Ə=<N 5. 
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The second approximation can be understood as neglecting the 
correlations between electrons located on different adsorption 
sites. Only the first approximation is more serious because it 
means we have neglected the possible correlated distribution of 
the adatoms on the substrate surface. In this way our 


Fig.2.The adatom local density of states for nickel 
substrate for V equal 2 eV,3 eV,5 eV and 6 eV 
for panels A,B,C and D,respectively.The first 
curve in each panel corresponds to coverage 6-20, 
second for @=0.5 and last for 6&-1,respectively. 


calculations can be comparable with those done, for example, 
within the single-site coherent potential approximation 
[35,37]. The second approximation is fully justified for our 
system, as it corresponds to neglecting the interatomic 
electron correlations. in Hubbard type Hamiltonian. The 


Electron Correlation Effects... 329 


numerical calculations have been performed for the electron 
charge of the hydrogen chemisorbed on the model substrates 
intended to describe nickel and titanium metals (for details 
seé Ref. (34]). For all coverages the value of q = <n 4> + <n, yb 
is significantly lower than results of Hartree-Fock 
calculations for the same set of parameters. The neutralization 
of the adatoms with increasing coverage was obtained in 


42202. - 01 
ENERGY ENERGY 


Fig.3.The same as in Fig.1 except for titanium substra- 
te for V equal 2 eV,4 eV,7 eV and 8 eV,respectively. 


agreement with the experimental results. The adatom local 
density of states depends on the coverage, too. For monolayer 
coverages the indirect interaction between adatoms through the 
substrate electrons is so large that can lead to two nearly 
separated energy bands. In Figs. 1-3 we have presi nted the 
adatom local density of states for three cases. Fig. . presents 


results for symmetric case, Fig. 2 for nickel substrate case 
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and Fig. 3 for titanium substrate case, respectively. In all 
cases a similiar situation appears, i.e., for large coverages 
(@ = 1 corresponds to one monolayer coverage) and for greater V 
one can observe appearance of a two electron subbands — effect 
very visible in Fig. 1. It can be understood as a result of the 
indirect interaction between adatoms through the substrate 
electrons which for large V and © is extremely strong. 


4. GENERALIZED HARTREE-FOCK SOLUTIONS 


The model Hamiltonians considered in preceeding sections 
have an unusual structure. They have the form of the Ising 
model with respect to ion variables Ny but with the operator 
-coefficients consisting of Fermi amplitudes [31]. It seems that 
the electronic properties described by correlators like SN a 
NN, se can be computed simultaneously within the same Green 
function "machinery" within equation of motion method. In the 
equation of motion method one obtains, in successive steps, the 
sequence of coupled equations for higher and higher order Green 
functions and the decoupling process usually applied in such 
Situations is an approximation, the meaning of which is not 
immediately transparent. For that reason, if the decouplings of 
the higher order Green functions into lower ones are made on a 
sufficiently high level of the hierarchy of equations or we 
take into considerations infinite series of the Green functions 
belonging to some classes, than we may believe that these- 
approximations can be justified. 

Here. we are going to calculate the Green function 
« Naso N35,» within the equation of motion method. Let 
us write equation of motion for this function: 


(E-€,)« Na Qe |Ng Vag? = Bc, < Na Ng) + 


(10) 
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where .]...» on the right-hand-side of the Green functions 

denotes IN 850”: Writing the equations for the Green functions 
appeared in (10) and for new Green functions in next equations 
we obtain: 


, | 
(E-e) & Na agl} VAES NaNe ne lise) an 


Ere 4 NN I T ba p C Na Na Ng? r 


ru & Na Na, Oa e la, e | eet 2 Vs i € Na Nec, Or, | He (12) 


(E- €) «NN rl >E Ver EN NN Dus, 0» 


(E- €,,) Naa, Na, |...» = Ap C Na Na, Nag Np) 

*U€N IN Neu gs Rag- |... * 22 Naat € NN Nee Oe |...) ada 

FE eae Na Ne... Ng; Qc |... "ap NATI Nap Ng) + 

ULN T Nepli; ai-e |» UNE n Ne gu... (15) 
- k t n= m 


Having in hand Eqs. (10-15) we can obtain exact equation for 


, * ) h - 
the Green function « Ngao 5525,» in terms of ionic correlators 


<N ANY (NS? and the Green function belonging to the specific 
a 


class of higher order Green functions, namely to the class of 
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where P on E= i D 


k 
Let us introduce some approximations in further treatment 


of this equation. First of all we approximate the  higher-order 
Ising-type correlators in the spirit of the Kirkwood 
superposition formula: 


QN. Na, +=- Nan? = CNa Na, --- CN, a Naa YL Na Na /4 AULA 


which was first used by Bethe [32] and Peierls [33] in the 
theory of lattices [38]. In this case the second term on the 
right hand side of Eq. (16) can be summed up exactly 
and reads as: 


4. iG (a-p) 4 
C Na Np) N a RE 4 (18) 


where P(E) is a Fourier transform of the Pap (P 


P j- Rep (E) Na Np > 


«p El (E-6) Na) 


It should be remarked that in our notation -vectors are 


(19) 


1 


two-dimensional vectors corresponding to the array of the 
adsorption centres in the adlayer and K-vectors are 
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three-dimensional. In a second step we neglect the last term 
of Eq.(16) containing the higher order Green functions and the 
third term decouple within an usual scheme. Note. that 
neglecting of the last term does not mean that we do not work 
with the infinite system of the equations for higher-and-higher 
order Green functions. This denotes simply that we removed from 

our considerations only special kind of higher order Green 


+ 


functions, namely, such as « NaNa coo NS 2a -o Nos. and all 
others present in Eqs. (10-15) up to infinite degree of 
hierarchy of this equations are present in our calculations. To 
this end, for the single-site Green function (a=/3) we obtain: 


UXN, "a.c 
€ Na due [Ot - 1; 27 | E -Ea ERO NE " (20) 
Y 


4 


- (1-8) A, (E) -8R; (8, £)] 


where we have used approximation (9). Here, ALE) and EE E) 
are as follows 


3 € Na, Na -6? 


(21) 


R (GE) =R (E) (4- SOT), on Qua) 


The Green function in the form represented in Eq. (20) 
generalizes the Hartree-Fock result in which ASQ(CE) and P'($,E) 
are replaced by A(E) and P(4,E) (38). 

As an another treatment of the system of equations (10-11) 
we take into consideration all the higher order Green functions 
which were omitted in our previous calculations. Unfortunately, 
we are not able to obtain an exact soliütuyons Of course Ror 
that reason we have to do some additional  *pproximitions. We 
decouple the higher order Green functions as follows (it is 
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generalization of usually used decouplings) : | 


<N Na, .-- Na,» 
CN Ns, Nes Qa gag | 2299 7 5 x | 
Ne Nan N 
{No Gags Ta, se | D = vum OE » (21) 
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Using Eq.21 we can rewrite two last parts of Eq.(16) in the 
form: 


6 
UC ne-e A $T Pigs SIE E (22) 
8 pez 4- Py (E) 


After straightforward calculations we obtain: 
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Using the Bragg-Williams approximation for 


calculation of 
P+ (E) and Bs (E) we have the final result: 


ew E 0) 
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where 
cial ry [E- €,-(-6) A (E)-6Pg (€)} ‘= 
ae s 6; HF(U =0) igh 
G 


pon E-&- E etae) - (1-8) A (E) -8 Pg (E) 


This formula stands for generalization of the. Hartree-Fock 
result for final coverages. All the higher-order Green 
functions appearing in the equation of motion for 
Na, ao Ngano? enter through two first terms in Eq. (16) 
Corresponding correlators are decoupled along the Kirkwood 
superposition formula.and ali other Green functions belonging 


to the class of functions of a type « NoMa: Na ĉa-o?a, -o Noaho”? 


are decoupled within the same method, i.e.. Kirkwood 
superposition formula combined with Hartree-Fock type 


decouplings. 


5. CONCLUSIONS AND REMARKS 


We have investigated the generalized Newns-Anderson 
Hamiltonian intended to describe the hydrogen atom or hydrogen 
overlayer chemisorbed on the metal surface. This Hamiltonian 
was obtained from more general one derived from microscopic 
considerations and'can be served as “ab initio" derivation of 
the chemisorption model Hamiltonian. We have retained only this 
terms which are important for our „specific, systems. In this 
case we have obtained the standard Newns-Anderson Hamiltonian 
supplemented by the additional terms describi.g the influence 
of the fractional occupation of the adatom inducec resonances 
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on the charge transfer between the adatom and substrate 
surface. The electron correlation effects have been considered 
beyond the Hartree-Fock approximation along the 
second-order self-energy matrix formalism of Ref.[16]. We have 
used model in which the nonorthogonality of the adatom and 
substrate wave functions is not taken into account, but in Ref. 
{4] has been shown that overlap effects nearly do not influence 
the adatom charge. On the other hand, the overlap effects may 
be essential for the calculation of the binding energy. For 
that reason we do not calculate the binding energy .The 
Hartree-Fock calculations performed for standard Hamiltonian 
predict the charge transfer to the adatom not compatible with 
experimental evidence. In our case results for generalized 
Hamiltonian cre much better. Even in the Hartree-Fock 
description calculated adatom electron charge is more realistic 
and after introduction of the correlation effects the results 
are omuoh better. In a next step we have considered a 
submonolayer coverages within another variant of a general 
Hamiltonian.Now, having in mind such complicated system we 
restricted ourselves to NA model with possibility of 
Simultaneous occupation of the adsorption centres by the 
adatoms.Using the equation of motion for appropriate Green 
functions and Bragg-Williams approximation for Ising type 
correlators we were able to take into account the stochastic 
arrangement of the adatoms within the overlayer.The nature of 
such approximation (in treatment of the overlayer) should be 
comparable with the coherent potential approximation but is 
much simpler in practical calculations.The variation of the 
electron charge localized on hydrogen adatom with increasing 
coverage is rather small „although ‘leads to progressive 
neutralization.This electron charge is too tlarge,of course, 
irrespective of the improved method along which the correlation 
effects were introduced.But this fact is rather result of the 
Simple model Hamiltonian (without additional term 

-IN V * + h.c. ) which was used. Taking into 


Ow O aanak a-o eo? ko 
consideration this term we should obtain improved results 


a 


similarly as in the one-adatom case [30].From the other hand, 
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these relatively small changing of the adatom electron charge 
vs. coverage are in good agreement with recently published 
results of the " ab initio " calculations [14]. Additionally, 
we ‘have also analyzed the adatom density of states for various 
coverages of randomly arranged adatoms.Generally one can 
say,that with increasing coverage a some sort of metallization 
of the overlayer can occur.In last section we derived the 
general expression for the one-particle Green function 
Naag Ngan,” for final coverages of the metal surfaces.We hope 
this general formula, Eq. (24), better describes the problem 
under cosideration as it takes into account all contributions 
although calculated in approximate manner.The analysis of the 
presented results calculated on the basis of the generalized 
Hamiltonian and the results given in Refs.[5,30,34] show that 
there is a satisfactory description of the chemisorption 
process. Thus,this Hamiltonian can serve as a model more 
efficient in describing the chemisorption process. 
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